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1 Introduction and main results 


The six-vertex model and related spin-i XXZ chain play a central role in the theory of exactly 
solved lattice models Q. Typically the six-vertex model is ‘solved’ by diagonalising the row- 
to-row transfer matrix with periodic boundary conditions. Several methods have evolved for 
doing this, including the co-ordinate Bethe ansatz Si, the algebraic Bethe ansatz |^, 
and the analytic ansatz ||5|. All of these methods rely heavily on the conservation of arrow 
flux from row to row of the lattice. 

In terms of the vertex weights (see figure 1) 


a = p sinh — v), b = p sinh i(A-|-u), c = p sinh A 

the transfer matrix eigenvalues on a strip of width N are given by 

A{v)q{v) = (/>(A — v)q{v + 2\') + 4>{X + v)q{v — 2X') 

where 

X' = A — ivr 
(f>{v) = /9^sinh^(iu) 

n 

<i{v) = 

k=l 

The Bethe ansatz equations follow from ( |1.2D as 

- Vj) ^ q{vj - 2A0 . ^ 

(j){X + Vj) q{vj + 2X')' ^ 

The integer n labels the sectors of the transfer matrix. 

Here we consider the same six-vertex model with anti-periodic boundary conditions. That 
such boundary conditions should preserve integrability is known through the existence of 
commuting transfer matrices [^J. However, the solution itself has not been found previously. 
In section 2 we solve the anti-periodic six-vertex model by the ‘commuting transfer matrices’ 
method |]I|. This approach has its origin in the solution of the more general 8-vertex model 
0, which like the present problem, no longer enjoys arrow conservation. We find the transfer 
matrix eigenvalues to be given by 


( 1 . 1 ) 

( 1 . 2 ) 

(1.3) 

(1.4) 

(1.5) 

( 1 . 6 ) 


A{v)q{v) = (/)(A — v)q{v 2X') — 4>{X -|- v)q{v — 2X') 


(1.7) 


where now 


N 

(l{v) = n sinhi(u - Vk). 

k=l 


( 1 . 8 ) 
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In this case the Bethe ansatz equations are 

4>{X — Vj) q{vj — 2X') 




(1.9) 


(/>(A + fj) q{vj + 2X'y 

In contrast with the periodic case the number of roots is fixed at N. 

In section 3 we use this solution to derive the interfacial tension s of the six-vertex model 
in the anti-ferroelectric regime. Defining x = e“^, our final result is 

' \2 


^-s/ksT _ 


2x5 


m=l 


1 -|- x‘ 


4m—2 


( 1 . 10 ) 


in agreement with the result obtained from the asymptotic degeneracy of the two largest 
eigenvalues ||l], 

With anti-periodic boundary conditions on the vertex model, the related XXZ Hamilto¬ 
nian is 

N 

i=i 

where and cr^ are the usual Pauli matrices, with boundary conditions 


•^TV-l-l — : 


y y z z 

^N+l ~ “‘^1) 


( 1 . 12 ) 


This boundary condition has appeared previously and is amongst the class of toroidal bound¬ 
ary conditions for which the operator content of the XXZ chain has been determined by 
finite-size studies Q. It is thus an integrable boundary condition, with the eigenvalues of the 
Hamiltonian following from (1.7) in the usual way Gt], with result 


E = NcoshX - V - 


^ 2 cosh sinh A 


(1.13) 


sinh ^Vj -|- sinh ^A 

We anticipate that the approach adopted here may also be successful in solving other 
models without arrow conservation. The solution given here can be extended, for example, 
to the spin-S generalisation of the six-vertex model/XWZ chain |10|. 


2 Exact solution 


To obtain the result (^) we adapt where appropriate the derivation of the periodic result 
Q (specifically, we refer the reader to Ch. 9 of Ref. |]l|). We depict a vertex and its 
corresponding Boltzmann weight graphically, as 


w{n,a\y,ij') = M- 
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where the bond ‘spins’ /i, a, (3, n' are each +1 if the corresponding arrow points up or to the 
right and —1 if the arrow points down or to the left. Thus in terms of the parametrisation 
O the nonzero vertex weights are 

n;(+,+!+,+) = = a 

'u;(+, —| —,+) = ui(—,+1,+, —) = 6 (2-1) 

w{+:-\+,-) = '“^(-,+1,-,+) = c 


The row-to-row transfer matrix T has elements 

01 02 

^a/3 = E • • • E - 

Ml MiV 


0N 


MiV 


• Miv+i 


( 2 . 2 ) 


^2 

«1 OL2 OLl\l 

where a = {ai,..., a^}, (3 = {/3i,..., Pn}, and the anti-periodic boundary condition is such 
that fJ,N+i = — 0 - 1 - Now consider an eigenvector y of the form 

y = gi® 92'^ ■ ■ ■ ® 9N (2.3) 

where ypap are two-dimensional vectors. From ( p.2| ) the product Ty can be written as 

{Ty)a = TV [Gi(ai)G2(a2) ■ • • G'Ar(a,v)5] (2.4) 

where Gi(±) are 2x2 matrices with elements 


0 


0 

The appearance of the spin reversal operator 

5 = 


gpp)- 


(2.5) 


0 1 
1 0 


( 2 . 6 ) 


in (2.4) is the key difference with the periodic case. However, it does not effect much of the 
working. Using (2.1) and (p.5|) we have 


Gi(+) — 


agi{+) 0 


Gi{-) = 


bgp-) cgi{+) 


cgi{-) bgi{+) ) \ 0 agi{-) 

In particular, there still exist the same 2x2 matrices Pi,..., Pn such that 

GPa) = 


(2.7) 


( 2 . 8 ) 
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where Pi and Hi are of the form 


V M-) *; Vo j 


(2.9) 


As for the periodic case, (|2.8| ) follows from the local ‘pair-propagation through a vertex’ 
property, i.e. the existence oi gi{a), gl{a),pi{a),pi+i{a) such that 


0P 

a\f3,p)gi{l3)pi+i{p) = gi{a)pi{p) 

(2.10) 

for a, p = ±1. The available parameters are Q 


9i{+) 

= 1, 5,(-)=r,e(^+’^)"«/2 


9i{+) 

= a, p'(-) = 

(2.11) 

Pii+) 

= 1, Pi{-) = P 


where cjj = ±1 and 


(2.12) 


However, pat+i needs to be different from the periodic case (where pat+i = Pi)- The anti¬ 
periodicity suggests that we require 


/WT)V jM-)\ 
V Piv+i(-) ) V p^(+) J 


(2.13) 


where h is some scalar. Since we already require Pi{+) = 1 and Pi{—) = Tj, we must have 

ri = — = —r and ttv+i = h = -. (2-14) 

h r 


In addition. 


^2 _ ^A((ti+...+(T jv) 


(2.15) 


To proceed further, we write Pi and Pn+i in full, 

A = { = 

Pl(-) ?l(-) 


Then 


hpi{-) qN+i{+) 
hpi{+) qN+i{-) 


0-1 

N-\-l 


SPi = 


1 

( qN+i{-) 

detPAT+i 

y -hpi{+) 

V l/h 

* 1 

0 - 

detPL 

detPiv+1 / 


(2.16) 


(2.17) 
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Putting the pieces together we then have 


(Ty)^ = TT[PiHiiai)---HN{aN)Pi^US 

= ■■■ 9N(.OiN) - h-^^—g'{{ai) ■■■ g%{aN). 


However, as for the periodic case, we have 


n, ^ ,gf{ai)detPi+i 
Qi {ai) = ah 

g[{ai)deiPi 


which follows from (2.7) to (2.9). Thus 




(2.18) 


(2.19) 


( 2 . 20 ) 


At this point it is more convenient to write 


y{v) = hi{v) 0 h 2 {v) 0 • • • 0 hjq{v) 


( 2 . 21 ) 


where we have dehned 


hi{v) = 


1 

n e 2 


The result ( 2.20| ) can then be more conveniently written as 

1 


T{v)y{v) = —ra^y{v + 2A') H— h^y{v — 2A'). 

r 


Also let 


N 

vtiv) = exp(-^^o-i) y{ 
^ i=l 


aj 


with r = ^exp(Y J2iLi o'i)- Then from ( 2.23 ) we have 

Tiv)yt{v) = ±(piX - v)yt{v + 2A') =f <('(A + v)y^{v - 2A'). 


( 2 . 22 ) 


(2.23) 


(2.24) 


(2.25) 


To proceed further, let Q^{v) be a matrix whose columns are a linear combination of y^ 

iN 


± 


with different choices of a {2^ altogether). It follows from ( 2.25 ) that 

T{v)Q'^{v) = ±(/>(A - v)Q^{v + 2A') =f <P{X + v)Q^{v - 2A')- 


(2.26) 


One can show that the transpose of the transfer matrix has the property T{—v) = ^T{v). 
With Q^{v) = ^Q^{—v) it follows from ( 2.26 ) that 


Qi{v)T{v) = ±4>{\ - v)Q^{v + 2A0 T </)(A + u)Q±(u - 2A'). 


(2.27) 
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Now let Qr{v) = Q^{v) and Ql{v) = Then we can show that the “commutation 

relations” 

Ql{u)Qr{v) = Ql{v)Qr{u) (2.28) 

hold for arbitrary u and v. This result follows if we can prove that {~'u)y^i{v) 

is a symmetric function of (w, u) for all choices of a, a'. Using (2.24), (2.21) and (2.22) this 
function reads 

\ 2 = 1 2=1 


3 = ^ 

X exp <{ i(A-u)o-j + ^(A + u)o-'- 


N 


i-1 


Y^{(yi+<T[) - ^{cri + <T[) 


. (2.29) 


_i=] i=l 

Now suppose that in cr, a' there are p pairs (cr^, cr'^) where = 0 with k = 1,... ,p. 

The terms in which involve these ai^ (in the prefactor and in the j = ik terms) are 
manifestly symmetric in {u,v). The remaining terms are exactly of the form (2.29) with 
N ^ N — p after relabelling of sites. We can thus restrict ourselves to the case where 
CTj = cr), i = 1,... ,N', for all N'. To prove this case we proceed inductively. From (2.28) we 
have 

N 


F = 

J- rrrr 


g-A((Tjv+...+o-j+i) gA((Tj_i + ...+o-i) 


j=i 


(2.30) 


Let us now denote Fq-o- = F]\f{ai,... ,aN)- By inspection, Fi(cri) and F 2 {ai,a 2 ) are sym¬ 
metric in {u,v). Suppose F]y- 2 {cri, ■ ■ ■, o'n- 2 ) is symmetric in {u,v) and furthermore that 
(^k + CTfc-i-i = 0 for some k. Then from (|2.30 ) we have F/v(cri,..., cr^, —cr^,..., cjtv) = 
Fjv- 2 (cri,..., (Tfc, (Tfc+i,..., crjv- 2 ) times a symmetric function of (w, v), which is therefore sym¬ 
metric in (u, v). This is true for all 1 < A: < — 1. The only case left to consider is therefore 

<7i = cr 2 = ... = cttv- But from (|2.30[) we have FAr((Ti, (T 2 ,... ,crAr_i,cri) = Fn- 2 {o' 2 , ■ ■ ■, o'at-i) 
times a symmetric function of {u, v), which is again symmetric. Thus by induction on N, the 
follows. 


assertion 

As in the periodic case, we assume that Qr{v) is invertible at some point v = vq and 
define 

Q{v) = Qr{v)Q^{vq) = QI^{vq)Ql{v). (2.31) 

Then from ( 2.27 ) and ( 2.2^) we obtain 


T{v)Q{v) = Q{v)T{v) = 4>{X - v)Q{v + 2A') - 0(A -|- v)Q{v - 2A') 


(2.32) 


‘Equivalent results are obtained using the other choice of sign. 
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and Q{u)Q{v) = Q{v)Q{u). This allows T{v),Q{v) and Q{v ± 2A') to be simnltaneously 


diagonalised, yielding the relation (1.7) for their eigenvalues. The precise functional form of 
the eigenvalue q{v) of Q{v), given in (1.8), follows from ( 2.32| ) by noting that T{v + 27ri) = 
—T[v) and considering the limits v —>■ ±oo. 


3 Interfacial tension 


In this section we derive the interfacial tension by solving the functional relation (1.7) and 
integrating over the band of largest eigenvalues of the transfer matrix |T^] . We consider the 
case where N, the number of columns in the lattice, is even. The partition function of the 
model is expressed in terms of the eigenvalues A(r;) of the row-to-row transfer matrix T{v) as 


z=i: iA(ioi 


M 


(3.1) 


where the sum is over all 2^ eigenvalues. 

The interfacial tension is defined as follows. Consider a single row of the lattice. For a 


system with periodic boundary conditions, in the A —> oo limit we see from (1.1) that the 
vertex weight c is much greater than the weights a and 6, so in this limit, the row can be in 
one of two possible anti-ferroelectrically ordered ground states. These are made up entirely 
of spins with Boltzmann weight c, and are related to one another by arrow-reversal. 

When we impose anti-periodic boundary conditions, this ground-state configuration is no 
longer consistent with N even. To ensure the anti-periodic boundary condition, vertices with 
Boltzmann weight c must occur an odd number of times in each row. Thus the lowest-energy 
configuration for the row in the A ^ oo limit will consist of — 1 vertices with weight c, and 
one vertex of either types a or b. This different vertex can occur anywhere in the row. 

As we add rows to form the lattice, the a or 6 vertex in each row forms a “seam” running 
approximately vertically down the lattice; it can jump from left to right but the mean direction 
is downwards]^ A typical lowest-energy configuration is shown in figure 2. The extra free 
energy due to this seam is called the interfacial tension. This will grow with the height M of 
the lattice, so we expect that for large N and M the partition function of the lattice will be 
of the form 

Z ~ exp [{-NMf - Ms)/kBT] (3.2) 


where / is the normal bulk free energy, and s is the interfacial tension. 


^ This is the analogue of the anti-ferromagnetic seam in the Ising model |l^. 
















We introduce the variables 


X = e 


-A 


z = e-^/^. 


(3.3) 


Expressing the Boltzmann weights in terms of z and x, from (O) the model is physical when 
z and X are real, and 2 ; lies in the interval 


X 


< . < .- 1 / 2 . 


(3.4) 


We consider A > 0 in order that the Boltzmann weights are non-negative, so we must have 
X <1. Let 


N 

Q{z) = “ ^3) 


(3.5) 


i=i 


where Zj = e j = 1,..., Ai, and 


V{z) = A{v){2zp-^)^ . 


(3.6) 


In terms of these variables, the functional relation (1.7) becomes 

Q{z)V{z) = (1 — z^x~^)^Q(—zx) — (1 — z'^x)^Q{—zx~^). 


(3.7) 


Both terms on the right hand side of (^) are polynomials in 2 ; of degree 3N, but the 
coefficients of 1 and z^^ vanish, so z~^V{z) is a polynomial in z of degree 2N — 2. We know 
how to solve equations of this form for both V(z) and Q{z) using Wiener-Hopf factorisations 
(see references 0, |] and ill ). 

We shall need some idea where the zeros of the polynomials Q(z) and V{z) lie in order 
to construct the Wiener-Hopf factorisations. From the anti-periodicity of T{v) we see that 
V{z) is an odd function of z, 

V{-z) = -V{z) (3.8) 


so its zeros and poles must occur in plus-minus pairs. To locate the zeros in the z-plane, we 
consider 2 ; to be a free variable, and vary the parameter x, in particular looking at the limit 
X ^ 0. 

We find the following; in the x —> 0 limit, iV — 2 of the N zeros of Q{z) lie on the unit 
circle, the other two lying at distances proportional to x^/^ and x“^/^. For V(z), there is the 
simple zero at the origin, and two zeros on the unit circle. The remaining 2N — 4 zeros of 
V(z) are divided into two sets, with iV — 2 of them that approach the origin and N — 2 that 
approach 00 as x —> 0. The N zeros of the two polynomials that lie on the unit circle are 
spaced evenly around the circle. 
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As X is increased, the zeros of Q{z) and z~^V(z) will all shift. We assume that the 
distribution of the zeros mentioned above does not change significantly as x increases. Thus 
the zeros that lie at the origin in the x —> 0 limit move out from the origin as x increases, 
but not so far out as the unit circle, and similarly for the zeros that lie at oo. Also, the zeros 
that lie on the unit circle are assumed to stay in some neighbourhood of the unit circle as 
X increases (we will show that these zeros remain exactly on the unit circle as x increases, 
which is what happens in the periodic boundary condition case). 

Bearing in mind the above comments, we write 

Q(z) = Qi(z)(z - a)(z -/3~^) (3.9) 

where Qi(z) is a polynomial of degree N — 2 whose zeros are 0(1) as x —> 0, and a^(3 = 
0(x^/^), so a lies inside the unit circle, (5 ^ outside. 

Also, let V{z) = z{z — ti){z — t 2 )A{z)B{z), where A{z) and B{z) are both polynomials 
of degree N — 2, the zeros of A{z) being all the zeros of V{z) that lie inside the unit circle, 
B{z) containing all those that lie outside, and ti and O are the zeros that lie on the unit 
circle. Since V{z) is an odd function, both A{z) and B{z) must be even functions of z, and 
we must have ti = —t 2 , so letting ti = —t 2 = t, we write 

V{z) = z{z^ — t^)A{z)B{z). (3.10) 


Draw the contours C+ and C_ in the complex z-plane, both oriented in the positive 
direction, with C_ outside the unit circle, C+ outside C_, and such that there are no zeros of 
Q{z) or V{z) on the boundary of or inside the annulus between C_ and C+. Then and 
all the zeros of B{z) lie outside C+ (see figure 3). 


Define r{z) as the quotient of the two terms in the RHS of the functional relation (3.7); 

Q{—zx~^){l — z'^x)^ 


r{z) = -- 


(3.11) 


Q{—zx){l — z‘^x~^)^ 

{r{z) has no zeros or poles on or between the curves C+ and C_). Then in the x —> 0 limit, 
we see that \r{z)\ ~ '^/z^, so when \z\ > 1, \r{z)\ < 1. Thus ln[l + r{z)] can be chosen to be 
single-valued and analytic when z lies in the annulus between C- and C+. We can therefore 
make a Wiener-Hopf factorisation of 1 -|- r{w) by defining the functions P+{z) and P-{z) as 

ie± " - z' - z 


lnP±{z) = ±^£ ln[l + r{z')] ^ 


(3.12) 


Then T+(z) is an analytic and non-zero (ANZ) function of z for z inside C+, and P-{z) is an 
ANZ function of 2 for 2 outside C_. As jz] ^ 00 , we note that P-{z) —> 1. When z is inside 
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the annulus between C_ and C_|_, we have, by Cauchy’s integral formula 


1 + r{z) = P+{z) P-{z) = 


V{z)Q{z) 


Q{—zx){l — z‘^x~^)^ 

We then define the functions V±(z); 

v+i^) = P+{z)Q{-zx)/{z - P~^) 

V-{z) = P-{z){l - z'^x~^)’^ / Qi(z){z - a) 


(3.13) 


(3.14) 

(3.15) 


where V+{z) is an ANZ function of z for z inside C+, V-(z) an ANZ function of 2 ; for 2 ; outside 
C-. We have split V{z) into two factors, l+( 2 ;) and V-(z), with 1 ^( 2 ;) = V+{z)V-{z) when 2 ; 
is between C+ and C_. 


Equating (|3.1C| ) with the expression 1 ^( 2 ;) = V+{z)V-(z) we have 

-z{z^-t^). 


P+{z) _ ^{z) ^,2 + 2 \ 


(3.16) 


B{z) V.{z) 

The LHS (RHS) is an ANZ function of 2 ; inside C_|_ (outside C_), which is bounded as | 2 :| —> 00 
and so the function must be a constant, ci say. Thus 


14(2;) = ciB{z) 

V-{z) = Ci^z{z^ — t^)A{z). 


(3.17) 

(3.18) 


When 1 2 ; I < 1, we proceed the same way. Draw the curves C(,_ and C'_, C'_^_ inside the unit 
circle, C'_ inside C(|_, and with a and all the zeros of A( 2 ;) inside C'_. 

In the limit x ^ 0 , |l/r(z)| ~ z^, so |l/r( 2 ;)| < 1. Thus ln[l + l/r( 2 ;)] can be chosen to 
be single-valued and analytic between and on C(|_ and C'_. We can then Wiener-Hopf factorise 
1 -|- 14 ( 2 ;) by defining the functions PLiz) and P'_{z) as 


1 


\nP'^{z) = ± — 


27 rz Jc' 


In 


1 + 


dz' 


r(2;')J z' — z' 


(3.19) 


where P\.{z) is ANZ inside C(,_, PL{z) is ANZ for 2 ; outside C'_. As \z\ 00 , PL{z) 1. 

When 2 ; is in the annulus between C, and C'_, Cauchy’s integral formula now implies 


1 + J_ = P' (,) p: U) = ___r(£)«£)_ 

r{z) Q{—zx~^){l — z‘^x)^ 


Define VL{z) and V'_{z) as follows: 


Vl{z) = P'^{z){l-z^x)^ I Qi{z){z -(5 ^) 
V^(z) = P!_iz)Q{-zx~^)/(z - a). 


(3.20) 


(3.21) 

(3.22) 
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We have now factorised V{z) into two factors, V\.{z) which is is ANZ for z inside and 
W (z) which is ANZ for z outside C'_. When z is in the annulus between and C'_, we have 
the equality V{z) = V!^{z)V!_{z). 

When z is inside this annulus, we equate (p.lO|) with V{z) = V^+(z)W (z) to get 


n(^) _ zAjz) 

B(z)(z2-t2) W(^) 


(3.23) 


where now the LHS (RHS) is an ANZ function of z for z inside C(|_ (outside €'_)■ Thus both 
sides of the equation are constant, C 2 say, and we have 


v:;(z) = C 2 {z^-t^)B{z) 
VLiz) = c^^zAiz). 


From equations (^1.17 ), ( 3.24 ) and ( 3.18 ), ( 3.25 ), we have the following 

n(z) = (c2/ci)14(z)(z2-t2) 

W(z) = ici/c 2 )VLiz)iz^ 


(3.24) 

(3.25) 


(3.26) 

(3.27) 


To evaluate the constant ci/c 2 , consider ( ^.271) in the limit z —> oo; we noted earlier that 
P_(z), PL{z) ^ 1 as z ^ oo, so from (^), (|3.15|) and (p.22D we deduce that 


Ci/C 2 = 1 . 


(3.28) 


We may use equations (|3.26|) and ( 3.27 ) to derive recurrence relations satisfied by Q{z), 
which we can solve explicitly in the N ^ oo limit. 

From equations (^.141) , (|3.21j) and (|3.26|) , we deduce the recurrence relation 


m Q(-zx} = (1 - - “X' -'’P 


(3.29) 


(z2 — t2) ^^(z) 

valid when z is inside C(|_. In the limit N ^ oo, the and functions —1, so we find 
that Q{z) is given by 


Q(z)= Q(o) n Y 


1 — z 2 x' 


2™4m—3 


N 


m=l 


(1 - z2t-2x^’^-2) (1 - za-^x2"^-2) (1 - Zf3x^^-^) 

— z2x‘^™'“^ J (1 — z2t“23.4m-4j _j_ 


(3.30) 


This still contains the parameters t, a and (3. From (3.29) in the N —> oo limit, setting z = 0 
we note that 


n 2 


<5(0) = —t ^ 


(3.31) 
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From equations (p-lSj) , (|3.22|) and (|3.27]) , we get the recurrence relation 


Q{z)Q{-zx ) = {1 - z X ) 


- a){z - P ^)P-{z) 


(3.32) 


{z"^ — P) P'-{z) 

which is valid for z outside C_. Taking the limit N ^ oo once more, so that the functions 
P-{z) and P'_{z) 1, we get 

N 


Qiz )=n 


1 — z X 


—2 _4m—3 


(1 - z-^Px^^-^) (1 - z-^ax^'^-p (1 - 


— Z 2^4m Ij (1 — Z (1 + Z (1 + Z ^(3 

(3.33) 

To derive an expression for l^(z) valid between C_|_ and C(_, using equation (3.27), we have 


V{z) = V+{z)VL{z){z‘^ - P) 

= Q{-zx)Q{-zx~^){z^ -P) / (z - a)(z -/3“^) 


(3.34) 


We use (|3.30 ) for Q{—zx) and ( 3. 331) for Q{ —zx ^), and substitute into equation (|3.34D . This 
produces a lengthy expression for V{z) involving the parameters a,/3 and t, which simplifies 
when one considers the oddness of V{z). The poles of V(z) must occur in pairs, and this is 
only possible if a and (3 are related by 


a/3 = —X. (3.35) 

Substituting this in, the infinite products involving a and /3 cancel, and we get, from 


and ( 3.34 ) 


A(u) = G{z/t) {p/2x)^ 


m=l 


where 


G(z) = ±ix^l‘^{z-z-^) n 


m=l 


- z^x^™-! 

1 - Z- 2 x^™-l\^ 

(3.36) 

- z 2 x"^™+l 

1 - Z- 2 x 4 ™+l j 

/ 1 — z^x"^™ 

l_Z-2x4m 

(3.37) 

i 1 — 

1 — 2;“23;4m-2 I 


This expression for the eigenvalue is still dependent on the parameter t, different values of t 
corresponding to different eigenvalues of the transfer matrix. All we know about t so far is 
that it is bounded as x ^ 0, and that it lies on the unit circle in the x ^ 0 limit. We shall 
now show that it in fact remains exactly on the unit circle as x increases. 

We substitute into the functional relation (3.7), using equations 


and 


to get 


an expression for the product Q{z)V{z) which is valid when z is in the annulus between C+ 
and C'_. Substituting into (3.7), the function on the right hand side is equal to zero when z 
is one of the N — 2 zeros of Qi(z), or when z = ±t. For the latter case, substituting z = t 
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and —t, and dividing the resulting equations, we arrive at the following relation between a, x, 
and t 

(3.38) 


which means that t must satisfy 


Wi)l" = ±1 


where 4>{t) is given by 


Ht)=t n 


1 — t X 


2™4m—1 


m=l 


1 _^-2^4m-3' 


1 - I _ 


(3.39) 


(3.40) 


This implies that t lies on the unit circle for all x, there being 2N possible choices for t. The 
partition function depends on t only via so there are only N distinct eigenvalues. The 
right hand side of ( |3.7| ) also vanishes when 2 : is a zero of Qi{z) so in the same way we show 
that the zeros of Qi(z) lie exactly on the unit circle for all x. As the zeros lie exactly on 
the unit circle, we may shift the curves C- and C'_^_ so that they just surround the unit circle. 
Hence our expressions for Q{z) are valid all the way up to the unit circle; ( 3.30|) is valid for 
l^l < 1, and (p.33 ) is valid for \z\ > 1. 

We now evaluate the partition function, as defined in (p.l|), in the large-lattice limit. 
When V is real, the eigenvalues ( 3.36 ) are complex, so as A' 00 , the partition function, a 
sum over the N eigenvalues defined by ( 3.3S| ), becomes an integral over all the allowed values 
of t, 

Z = ^ p{t) [A{v)]^ dt (3-41) 

where the integral is taken around the unit circle, and p{t) is some distribution function. 


independent of N and M. Substituting (3.34) into (3.41) then gives an expression for Z. 
(The number of rows M is even to ensure periodic boundary conditions vertically, and so the 
± sign in ( p.37 ) is irrelevant.) 

The eigenvalue ( p.3(: ) contains two distinct types of factors; those that are powers of N, 
and those that are not. The terms that increase exponentially with N contribute to the bulk 
part of the partition function, the free energy per site in the thermodynamic limit. This 
factor is also independent of t, and can be taken out of the integral (p.41). The integral is 


then independent of A, so we have, from (3.2) 

00 

g-//fcsT ^ (^/ 2 x) f[ 


1 — ^ 


1 - z-‘^x^^-^' 


771=1 


^ _^2^4m+l 2^4m+l 


(3.42) 


for the free energy per site in the thermodynamic limit. This result agrees with the result for 
periodic boundary conditions (equations (8.9.9) and (8.9.10) of Ref. 0)- 
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From equation (3.2), the other factors in (3.34) make up the interfacial tension, given by 


^-MsIksT ^ ^ [G{z/t)]^ dt. (3.43) 

For M sufficiently large, we may evaluate this integral using saddle-point integration. The 
integral is given by the value of the integrand at its saddle point, together with some multi¬ 
plicative factors that we can disregard as M —> oo. The function G satishes the relation 

G{z)=Gi-l/z) (3.44) 

which implies that the function has a saddle point when z = ±i. Hence the integrand in 
(H) is maximised when 

i — ^saddle — iiz. (3.45) 

As z is arbitrary, this point may lie off the unit circle; it will however lie inside the annulus 
between C+ and C'_ because of the restriction (^^), and so we will be able to deform the 
contour to pass through this saddle point. Hence we arrive at the final result 

e-A-’' = 2,d/=n ■ P.46) 

m=l \ / 
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a a b b c c 

Figure 1: Standard vertex configurations and corresponding weights. 



Figure 2: A typical lowest-energy state of the system with N even and anti-periodic boundary 
conditions. The dotted line indicates the interface dividing the lattice into two domains, each 
of which is an anti-ferroelectrically ordered ground state. 
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